Attractive Hubbard model and single-particle pseudogap due to classical pairing 

fluctuations in two dimensions 
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Abstract: It is shown that in the two-dimensional at- 
tractive Hubbard model, the mean-field phase transition is 
replaced by a renormalized classical regime of fluctuations 
where a pseudogap opens up in the single-particle spectral 
weight. It is argued that this pseudogap and precursors of 
the ordered state quasiparticles can occur only in strongly 
anisotropic quasi two-dimensional materials. This precursor 
phenomenon differs from preformed local pairs. Furthermore, 
while critical antiferromagnetic fluctuations would also lead 
to a pseudogap in the repulsive model, there are some impor- 
tant differences between the superconducting and magnetic 
pseudogap. 

For several years, experimental evidence about the ex- 
istence of gap-like features in the normal phase of under- 
doped high-temperature superconductors has been accu- 
mulating. In the last year or so, photoemission experi- 
ments have convincingly shown that fermionic excitations 
do develop a pseudogap on the Fermi sur face §. This 
pseudogap, that forms well above the superconducting 
transition temperature, suggests a normal-state precur- 
sor of the superconducting gap. It has a similar angular 
dependence and basically the same magnitude. Several 
explanations of this gap have been put forward. They 
range from strong coupling effects, to magnetic effects, 
to preformed local pair effects 0|. In this paper, we 
present a Quantum Monte Carlo study of the attractive 
Hubbard model to show how precursors of the super- 
conducting quasiparticles may occur in the normal state, 
creating pseudogap features in the single-particle weight 
/ (a;) A (k,w) measured in Angular Resolved Photoemis- 
sion Experiments (ARPES). Pseudogap features caused 
by superconducting fluctuations have a very long history, 
dating back to film studies in the 1970's. || Such fluc- 
tuations have been proposed again more recently under 
the name of phase fluctuations Q as a cause of the pseu- 
dogap. The pseudogap features discussed before were in 
the total density of states. Conditions for a pseudogap in 
the momentum-dependent / (w) A (k,u/) are more strin- 
gent. These conditions have been studied in the magnetic 
context 0] and analogies with the superconducting case 
have also been pointed out. M 



In this paper, we show that classical thermal fluctua- 
tions are sufficient to obtain precursors of superconduct- 
ing quasiparticles at arbitrarily small coupling as long as 
a) the system is quasi two-dimensional b) the supercon- 
ducting correlation length £ grows faster with decreasing 
temperature than the single-particle thermal de Broglie 
wavelength £th = vp/T. This last condition is always 
satisfied in the vicinity of Kosterlitz-Thouless or mean- 
field transitions, although it is more stringent than the 
condition £ > a (a is the lattice constant) assumed by 
several authors. The Physics of this effect is very differ- 
ent from that of the so-called preformed local pairs that 
destroy the Fermi surface in any dimension but only for 
sufficiently strong coupling. 

The attractive Hubbard model, that leads to an 
s— wave state, is used here because it is easier to sim- 
ulate, allowing us to check the sufficient conditions ob- 
tained analytically for the appearance of precursors of the 
superconducting quasiparticles in the normal state. All 
simulations arc done for U = —At and units are chosen 
such that t = 1, a = 1, fcs = l,h = 1. 

At half-filling, where the chemical potential /i vanishes, 



maps 



the canonical transformation Cjj — > ( — 1) 
the attractive model onto the repulsive one. The q = 0, 
s— wave superconducting fluctuations and the q = (jr, w) 
charge fluctuations are mapped onto the three antifer- 
romagnetic spin components of the repulsive model and 
hence they are degenerate. Because of this degeneracy, 
the order parameter at half-filling has SO (3) symme- 
try and hence, by the Mermin- Wagner theorem, in two 
dimensions the phase transition cannot be at finite tem- 
perature. At half-filling the transposition of the results 
for the repulsive case shows that the pair structure fac- 
tor 5 A = ((At A) +<AAt» with A = -j=Y.ti^i 
becomes size-dependent and saturates at low tempera- 
ture. This saturation signals a superconducting ground 
state. Instead of a finite temperature phase transition, 
the sudden rise of 5a as temperature decreases indicates 
a crossover to a renormalized classical regime where the 
correlation length £ increases exponentially. In an in- 
finite system, a pseudogap opens up all the way from 
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the crossover temperature Tx to zero temperature. The 
opening of this pseudogap is signaled by a sharp de- 
crease with temperature of the following measure of the 
single-particle spectral weight near the Fermi energy: 

Away from half-filling, the SO (3) symmetry is bro- 
ken down to SO (2) by the finite chemical potential. 
The charge fluctuations are suppressed compared with 
the superconducting ones. The size-dependent data can 
be scaled like at half-filling with an exponentially grow- 
ing correlation length, showing that, for the accessi- 
ble temperature range, this is still an appropriate scal- 
ing. It is noteworthy that at the crossover temperature 
Tx the charge fluctuations are still comparable to the 
superconducting ones. Despite the fact that a finite- 
temperature Berezinski-Kosterlitz-Thouless (KT) tran- 
sition is possible with an SO (2) symmetry, the rapid 
rise of superconducting correlations apparent from Fig. la 
is a crossover to a renormalized classical regime, not a 
finite-temperature KT transition. |8) The fact that this 
is a crossover is demonstrated by Fig. lb which shows 
the ratio of the classical (zero Matsubara frequency) to 
quantum contribution to the pair structure factor. We 
see that the classical contribution decreases as tempera- 
ture decreases (from T = 2toT=l/2in the inset of 
Fig. lb) and then it starts to rise again at lower temper- 
ature, signaling the onset of the renormalized classical 
regime at Tx (n). The downturn at lower temperature 
indicates that as temperature becomes smaller than the 
minimal energy of the pairing fluctuations in the given fi- 
nite system, the finite system returns back to a quantum 
regime and eventually S&(L, q — 0) saturates to its zero- 
temperature value in the superconducting ground state. 



It is known from classical-spin simulations that the KT 
transition is seen in finite-size systems only for much 
larger lattice sizes than those accessible here, namely 
when at least several vortices can fit inside. 

In the renormalized-classical regime that exists be- 
tween Tx and Tkt one should observe features in the 
single-particle spectral weight that are precursors of the 
superconducting ground state quasiparticles. This is sug- 
gested by the drop in z (T) in Fig.lc that coincides with 
Tx- To show the precursors more explicitly, consider in 
Fig. 2 the single-particle spectral weight A(k — (0, ir) ;uj) 
obtained from analytical continuation of Monte Carlo 
data using Maximum-entropy techniques. Temperature 
decreases from T = 1/3 to 1/5 from the top to the bot- 
tom curve. These temperatures correspond respectively 
to the points marked A,B,C, in Fig.lc. Clearly, the 
Fermi liquid quasiparticle that exists at high tempera- 
ture, T = 1/3, starts to split into two superconducting- 
like quasiparticles around T = 1/4, opening a pseudogap 
at the Fermi level. The splitting at T = 1/4 occurs in 
a size-independent regime (6x6 near 8x8) while at 
T = 1/5 we start to be in the size-dependent regime of 
Fig. la. Note that k = (0, tt) is as close as possible to 
the Fermi surface. The pseudogap should be larger for a 
point strictly on the Fermi surface. 

The above results may be understood by analytical 
methods that have been validated in detail in the half- 
filled case. H Assuming some effective coupling con- 
stant g between quasiparticles and pairing fluctuations, 
one can obtain for the classical contribution to the spec- 
tral weight, (T > \ui + e_k| > «f£ _1 ), 

am*- --^0 a) 
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As long as £ increases exponentially, it is clear that 
A 2 oc gT In £ will be finite so that the spectral weight 
will have a maximum at frequency given by a BCS-like 



dispersion or 



-2 



A . Clearly the above effects can 



exist even in the weak coupling limit. This has to be 
contrasted with the case of preformed local pairs that 
appear when the coupling strength becomes of the or- 
der of the bandwidth. The latter is, certainly, not the 
case in High-T c materials where the gap is about SOmeV 
while the bandwidth is of the order W « leV. Also, 
preformed local pairs can occur for U > W in arbitrary 
dimensions. Here however, phase space considerations 
for the integral defining the renormalized classical con- 
tribution show that in three dimensions one does not ob- 
tain the singular form Tln£/ {u> + e_k) for ^2' that is 
necessary to obtain finite-frequency precursor quasipar- 
ticles. Nevertheless the imaginary part of £ in three di- 
mensions is weakly (logarithmically) divergent close to T c 
which leads to a pseudogap in the total density of states 
N(ui) = A(k, oj) . Note that approximations that 

include self-consistency in the Green's functions without 
the corresponding vertex corrections (FLEX like approxi- 
mations) fail to reproduce the precursor of superconduct- 
ing bands in the k— resolved spectral function A(k, ui) 

Physically, precursor effects exist because in the 
Kosterlitz-Thouless picture the magnitude of the order 
parameter is locally non-zero starting below a crossover 
temperature Tx that is larger than the transition temper- 
ature T KT . It is only the phase that is globally decorre- 
lated above T^t- Another way to understand the precur- 
sor effects is that the superfluid density and the gap are 
finite as T — > T and, hence, a two peak structure in 
A(kp ,u>) exists even as the phase transition point is ap- 
proached from the low-temperature side. This two peak 
structure should not immediately disappear when one in- 
creases the temperature slightly above Tkt- In quasi- two 
dimensional systems these high energy precursors should 
persist. By contrast, in an isotropic three-dimensional 
system, the gap vanishes at the transition point and there 
is no finite frequency precursors even though there may 
be a low-frequency depletion of the total density of states. 

It is tempting to speculate that the physical origin of 
the decrease of T c in the underdoped materials is sim- 
ilar to what happens in the attractive Hubbard model 
due to the fact that Tkt must go to zero when the num- 
ber of components of the order parameter is larger than 
N > 2 ( 50(3) in the attractive model). However, at 
present the size of the fluctuating regime of the phase 
transition in high-T c materials is controversial. Recent 
infrared experiments H] suggesting a pseudogap even in 
the overdoped regime are consistent with the mechanism 
discussed here. Clearly however, present ARPES exper- 
iments M show large backgrounds that are unexplained 
by the above approach. Finally, let us contrast precursor 



effects due to magnetic and superconducting fluctuations. 
Our criterion for the appearance of a pseudogap shows 
that even for the s-wave interaction the superconducting 
pseudogap should open up from the zone edge (smaller 
vf) as temperature decreases until the whole Fermi sur- 
face is gapped. This is, obviously, even more so for the 
d— wave like interaction, when one has nodes along the 
diagonal direction. In the case of an antifcrromagnetic 
pseudogap, the region near the hot spots is gapped, but 
a whole segment of the Fermi surface near the diagonal 
direction remains gapless even when one approaches a 
phase transition. 
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